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' Possible phase transitions between incompressible quantum Hall states and compressible three- 

, dimensional states are discussed for infinite-layer electron systems in strong magnetic field. By 

variational Monte Carlo calculation, relative stability of some trial states is studied. If the inter- 
layer distance is large enough, the Laughlin state is stabilized for the Landau- level filling v = 1/3 of 
each layer. When the inter-layer distance is comparable to the magnetic length, the Laughlin state 
becomes unstable against a different fractional quantum Hall state with inter-layer correlation, 
and/or against a three-dimensional compressible state. It is discussed how the quantum phase 
transition between them can be controlled in actual systems. 
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Although the fractional quantum Hall effect is originally a two-dimensional phenomenon [Q-Q , multi-layer systems 
have been studied as natural development Multi-layer systems show some unusual effects because of new 

additional degrees of freedom, namely the inter-layer Coulomb interaction and the inter-layer electron tunneling. For 
example, in the case of the simplest double-layer system, the " even denominator" fractional quantum Hall effect with 
the Landau level filling v = 1/2 =1/4 + 1/4 has been observed 0], and is interpreted in terms of Halperin's ^331 
q 1 state ||. 

One may then naturally ask what happens in infinite multi-layer systems. Experimentally, the integer quantum Hall 
effect is observed in a 30-layer GaAs/AlGaAs superlattice system jlOj and also in organic conductors (TMTSF) 2 X 
and BEDT-TTF salts Motivated by these works we study the ground state of the infinite-layer system 
■^j- • under a strong magnetic field and discuss quantum phase transitions between incompressible and compressible states. 

We can have the following naive expectation in advance of quantitative analysis: fn the limit of infinite inter-layer 
distance d s , layers are completely independent of each other and the Laughlin state is realized for appropriate fillings, 
e.g. v = 1/q where q is an odd integer. However for d s ~ Ib where Ib is the magnetic length, other fractional 
quantum Hall states can be realized because of the inter-layer Coulomb interaction. Furthermore three-dimensional 
states are also realizable with increasing electron tunnelings between layers. The fractional quantum Hall states are 
incompressible, but the last three-dimensional state is compressible because the Bloch states are formed perpendicular 
to layers and the Fermi surface should appear. Therefore a transition between incompressible and compressible states 
I is expected at a certain value of d s , provided that d s is a variable parameter. 

Based on the above expectation we introduce three typical trial states for the ground state with the lowest Landau 
level filling v = 1/3 and compute energies of them in order to find which one is the ground state. Under a strong 
magnetic field we restrict the single- particle states within the lowest Landau level, and consider only the Coulomb 
energy and tunneling energy. We calculate the Coulomb energy by the variational Monte Carlo (VMC) method for 
the two-dimensional trial states Jl^Jr?|] , while the Hartree-Fock approximation is used for the three-dimensional trial 
I> ■ state. 

C3 1 II. MODEL 

We consider a system with the area S of each layer and the height L z as shown in Fig. [j]. The total number N z 
of layers is given by N z = L z /d s . Then, with the total number N of electrons in the system, we define the electron 
number per layer N±_ by N± = N/N z . Taking into account the layer structure of the system, we assume the positive 
charge distribution 

N * N 

n fe (r)=^-^(z-Zd s ). (2.1) 
1=1 
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This positive charge density corresponds to the uniform distribution in each layer and the discrete distribution along 
the layer. Under a magnetic field B = Be z the Landau level degeneracy mo is given by mn = S/2xl B where 

1/2 

Ib = (ch/eB) is the magnetic length. The average filling factor of each layer v is defined as v = N±/m£>. Under 
the strong magnetic field with v < 1 we may deal with only the lowest Landau level (LLL) explicitly and neglect 
higher levels. We will discuss possible ground states in the subspace of the LLL in the following. 
Now we write down the Hamiltonian of the system with N — > oo as 

N e* 2 N f e* 2 1 f e* 2 

H = J2 H *i + E 17- —\ E / dr ^— H n " (r) + 2 / drdr 'lr~~ FT 6 (r) m (r ' } ' (2 - 2) 

i=l i<j 1 2 i=l J 1 1 1 J 1 1 

Here the one-body part is given by 

= (px* + ^ A ^ ( r *)) 2 + ^Prf + V ° ^) > ( 2 - 3 ) 

where m*j_ and m| are the effective masses in the layer and along the z-direction respectively. The squared effective 

charge e* 2 is defined by use of the dielectric constant e as e* 2 = e 2 /e. Vq (z) is the binding potential of electrons to 
layers (like a multi square well potential). We choose two types of configuration as follows: one is the circular shape 
and the other is the rectangular shape. These shapes are chosen for convenience of calculations. 



A. Circular Geometry 

Let us consider a pile of infinitely large layers (S — oo) with the circular symmetry and take the symmetric gauge. 
For the (x, y)-space this geometry is the same as the one that Laughlin chose [g]. We adopt for this geometry the 
periodic boundary condition in the z-direction. This geometry is convenient for description of two-dimensional states 
and also convenient for the VMC calculation If we assume that electrons do not transfer between layers, then 
the one-particle state is written as 



i) lm (r) = (2 m+1 nmH 



2m+2\- 1 / 2 



fexp 



'4/1 



0o {z - ld s ) , 



(2.4) 



where £ = x — iy is the complex coordinate in the (x, y)-space and (f>o (z — ld s ) defined by (2.14) is the normalized 
wave function localized at the layer I. Since the uniform positive charge is inconvenient for finite particles which we 
actually study numerically, we take an alternative positive charge density hb (r) which has the same distribution as 
that of electrons. Namely with P ({r^}) = ^({""i})! 2 being the probability distribution function of electrons, we 
require 



n b (r) 



(2.5) 



This positive charge distribution tends to the uniform one in the infinite-size limit. 

We now discuss the Coulomb potential with periodic boundary condition along the z-axis. Taking account of mirror 
images we obtain the potential as 



n— — oo 



[(r_L -r'J 2 + (z-z' + nL z y 



1/2- 



(2.6) 



The summation over n is divergent logarithmically. This divergence should be canceled by the attractive part in the 
final result. In this paper we replace it by the following one: 



V (r - r') 



,*2 



(r-L - ^ 



^ sin 



1/2- 



(2.7) 



This potential has the right periodicity along the z-axis and becomes equivalent to the Coulomb potential in the limit 
of L z — > oo. The Hamiltonian with finite N is represented as follows: 
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N N 

h = J2 h ^+Y, v ^ - r j) - E / drV ( r * - r ) hb w + 2 / drdr ' v (r ~ r,) hb (r) " 6 (r,) 

• n a- jl X 1 « *■' 



(2.8) 



i — 1 2<J 

The expectation value is given by 



(H)=Nn+j2 [ ( n dr v ^ - r ->) p » r ^) 

i<j \i=l / 
t,j \t=l / 



where 57 is a constant representing the single-particle energy. It is given by 

Sl=\u c + j dz<& (z) + Vb (*)J 0o (z) , 

where we have assumed that the tunneling along the z-direction is suppressed. 



(2.9) 



(2.10) 



B. Rectangular Geometry 

Next we consider a pile of rectangular (S = L x x L y ) layers. We impose the periodic boundary condition in each 
(x, y and z) direction and choose the Landau gauge. We obtain mu — L x L y /2irlg. In the case of three-dimensional 
states, this geometry is convenient. When there are tunnelings, the one electron state is written as 

cf> jk (r) = <f, ±j (r ± ) 0|| fc (z) , (2.11) 

where j (= 1,2, ■••,m£>) indicates the degree of freedom in the LLL, and k = 2m/ L z (n — —N z /2,—N z /2 + 
1, • • • , N x /2 — 1) is the wave number along the z-direction. The (rjj-dependent part <f>±j (rj_) is the Landau wave 
function [HI: 



L y lB\/TT 



1/2 



E eX P 

n= — oo 



(2.12) 



where Kj = 2"Kj/L y and = l 2 B Kj. For the z-dependent part <f)^ (z) we make the Bloch function by use of 0o as 

oo N, 



1 



IN, 



J2 ^e 4fcnd ^o(^«4~n z i,) 



n z — — oo n— 1 

where we take the following localized function: 



00 (*) 



1 



2 d 2 n 



1/4 



(2.13) 



(2.14) 



with e being positive infinitesimal. Under the tight-binding approximation tfijk (r) becomes the eigenstate of the 
one-particle Hamiltonian H^. We introduce the creation and annihilation operators, <zt fe and ajk, corresponding to 
cf>jk (r). Then the total Hamiltonian of this system is given in second quantized form by 

N e* 2 d 2 



" = EE ^ 



3 fe 



2F 7T 



2C(2) UW 



q#0 



p (q) p (-q) - p (o) ex P 



;2 2 2 j2 2 



(2.15) 



Here Ef, — — 2t cos A;d s is the eigenvalue of the one-electron state with the inter-layer transfer —t, ( (n) is the 
Riemann's zeta function, v (q) = Aire* 2 jq 2 is the Fourier coefficient of the Coulomb potential, and p (q) is the density 
operator defined by 



p(q)= / drip^ (r) tp (r) e" 



iq-r 



(2.16) 



where "0 ( r ) = Sifc a jk4'jk (r) is the electron field operator projected onto the LLL. 
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III. TRIAL STATES 



We introduce three trial states for the grou nd states with v = 1/3 of the Hamiltonians (2.8) and (2.15) . A s 
mentioned in the previous section, we take ( |2.8| ) for states without the tunneling (two-dimensional states) and ( 2.15| ) 
for states with the tunneling (three-dime nsio nal states). Here we remark on the index of position coordinates 

■ , N) of electrons 



position coordinates in the Hamiltonian (2 
to two-dimensional states 



The 

However when we refer 



have the index i (= 1,2 
we make the position coordinates have the layer index additionally because electrons 
belonging to different layers are distinguishable. Hence the position coordinates have the two indices I and p where I 
(= 0, 1, • • • , N z — 1) is the layer index and p (= 1, 2, • • • , N±) is the particle index in a layer. Indices {i} and {l,p} are 
related to each other by i = N±l + p. 

The first trial state is a two-dimensional one without the tunneling represented by 



N z 

*-i/3 =n[^=v3 ({&»] <!>({{%>}}), 



(3.1) 



(=1 



where \ an d $ describe the motion in the (cc, y)-space and along the z-direction respectively. They are given by 



x£=i/3 ({&!>}) = II (^p ~ &p') 3 exp 
p<p' 



4/ 



1 N± 



p=l 



(3.2) 



*({{^}})=nn^(^- w )- ( 3 - 3 ) 
i=i p=i 

Here ({{zi p }}) denotes (zn, ■ ■ ■ , Zi N± ; ■ ■ ■ ; z Nz i, ■ ■ ■ , z NzNj _) and ({& p }) does (£a, £; 2 , • • • , €w±)- We call the state 
( |3.l|) the " multi-Laughlin state". It is the direct product of the Laughlin state with v = 1/3 || in each layer and is 
expected to be stable for a large inter-layer distance: d s ^> Ib- 

The second trial state is also a two-dimensional one without the tunneling. It is represented by 

*£=V3 = xS^({{^}}), (3.4) 

where 

^s=i/3 = n n ~ n ^ ~ ^+1^) ex p 

/ — 1 P<f)' P-.p' 

We call this state the " (111) state" J(|. This state is expected to be stable for d s ^ Ib where the inter-layer Coulomb 
interaction becomes comparable to the intra-layer one. The numbers (111) indicate the exponents of the Jastrow 
factor which controls the strength of the electron correlation. The first " 1" is the exponent of the Jastrow factor 
between a given layer and the nearest lower layer, the second " 1" is that within the layer, and the last " 1" is that 
between the layer and the nearest upper layer. Generally we can define (A/iA) state for our system. For example (131) 
state is possible for v = 1/5, and (212) state is also possible for the same v. Physically (212) state should be unstable 
because the inter-layer correlation becomes larger than the intra-layer one in this state. In the case of v = 1/3 only 
the (111) state is allowed within the subspace of Halperin-type states. 

The third state is expected to be stable at d <J Ib when electrons transfer between layers. We call such a state the 
"itinerant state" and define it as follows. With the tunneling, Bloch waves are formed along the z-direction. Then 
the simplest candidate for the ground state is the one with the smallest kinetic (tunneling) energy. Namely the cosine 
band in the z-direction is filled up to the Fermi level |l9| , with the filling of the Landau level for each wave number 
being unity below the Fermi level and zero above it. We can write it by use of electron creation operators as 

rriD 

| itinerant) = JJ JJ o} fc |0), (3.6) 

3=1 \k\<k F 

where fcp is the Fermi wave number given by kp — vir/d s . Although we assume v = 1/3 in this paper, the itinerant 
state can be defined for all v G [0, 1]. 



1 N± 



(3.5) 
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IV. GROUND-STATE ENERGIES 



A. Multi-Laughlin State and (111) State 

We use the VMC method to calculate the expectation value of the Coulomb energy numerically for the multi- 
Laughlin state and the (111) state. Technical details are explained in Appendix. In the case of the multi-Laughlin 
state, (H) are reduced to the sum of the two-dimensional result simply because all layers are independent of each 
other, and the inter-layer Coulomb energy is exactly canceled. We obtain 

Ne L = (H) L Nn = 7V ± ]T / ( Y[dr xlp ] _f^_| x ^ 1/3 ({fc p }) | 2 



^e / (n^p) 

p<p' \ p / 

>-E/(n^ 

p,p' \ p 



e 



*2 



-IP 1 ±.lp 



j—Axt= 1/3 Hx^r/3 ({fU) I'- 



On the other hand such cancellations do not occur for the (111) state, and the calculation becomes more complicated 
because of the inter-layer correlation. The energy is given by 

Ne (in) = (H) {lu) -Nn 

ee+ee) / fnn*O v(r *- r ' v)| *^/» ({{r *>> )|a 

K l=l> p<p> l<l> p.p'j \ I p J 

4ee / (nn^4>W p - r ^^ (42) 

p,p' \ J n / 



Figure || shows the size dependence of the Coulomb energy of the multi-Laughlin state. The extrapolated value is 
in good agreement with that in the literature (l7). The dependence of the Coulomb energy for the (111) state on the 
inter-layer distance is shown in Fig.||. The calculation was performed for three cases of finite size systems: N±_ — 20, 
30 and 42 with N z — 5. We find that the energy decreases with decreasing inter-layer distance. The dependence on 
N±_ is small in the region of small inter-layer distance. 



B. Itinerant State 



We calculate the energy of the itinerant state for the Hamiltonian ( 2.15| ) with use of the Hartree-Fock approximation. 
Since we take the homogeneous trial state, the mean field is given by 

A (jk; j' k') — (oijkO'j'k') — (itinerant|aj fc aj'fe' |itinerant) 

= 6 jf 6 kk ,6(k F -\k\). (4.3) 



The mean field for the interaction part of the Hamiltonian (2.15) consists of the direct and the exchange interactions. 
Of these the energy due to the direct term is written as 

1 iV 2 

57? E v ^) (p(q))(p(-q)> = w^2 v (q) s <ix,gK,h 



2V S-f v"^v-x„u-v ~w, w 

q/0 q#0 



-2C(2), (4.4) 



2V TT 

where we used the relation (p(q)) = N5q ±tg 5 q , t h valid in the macroscopic limit. Here we have g = 

^j^mi>n x , jZ-mjjny^J and h = ^f-n z with n x ,n yi n z integers. The direct term cancels exactly with the Coulomb 

energy of the positive charge because of homogeneity of the itinerant state within each layer. Therefore the Hartree- 
Fock Hamiltonian is given by 
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n HF = J2J:E k a] k a jk 

3 k 

~ 2V E E a )i ki a ii fe i 9 ( kF ~ I fc2 1 ) E v ^ M tii k i > h k 2 ; q) M O2 ^2 , ii fci ; -q) , 



(4.5) 



J1J2 k u k 2 q/0 

where the matrix element for the exchange mean-field is given by 

M (#,/*'; q) = /" dr*J fc (r) (r) e*»" 

00 

= / / ^fc' fc-0 z +/!<5- VI L » eX P 

^ > n ' lT " j— j '+m D n, — -2%-q y r 

h n— — 00 

The energy of the itinerant state is obtained as 
E = m D ^2 E k 

\k\<k F 



l 2 B q\ e 2 d 2 s q 2 



iq x Xj 



(4.6) 



E«(«)EE fl ^- 6 (** - 1* - «* + ^i) ex p 

q#0 



fe ft. 



(4.7) 



The region of summation over k in the above formula is shown in Fig. 0. In order to derive the sum S we analyze 
the cases as follows: (1) S = for q z - h > 2k F , (2) S = ^{k F -{q z - h - k F )) for < q z - h < 2k F , (3) 
S = i^(q z — h + k F — (—kp)) for —2kp < q z — h < 0, and (4) S = for q z — h < —kp. As a result we summarize the 
four cases into a single expression given by 



S = ^2k F [ I 
2tt 



h\ 



2k f 



'(2fc, 



h\) 



(4.8) 



Therefore the energy per particle which contains the kinetic term is given in the thermodynamic limit by 



e HF = a 



2t 



2(2*)° ^Jo+ 

H F 

Q! 1 Skin H~ S-int ' 



27rq±dqj 



h+2kF , \q,-h\ 
dq z v {q) 1 

h-2k F V 



2k f 



exp 



I 2 

2 q± 



(4.9) 



At this stage we let e — > 0. It can be shown that the contribution of large q z (~ l/ed s ) to the integration is negligible. 
To see this we observe that v (q) = 1/ (<j^ + 9 2 ) becomes (e 2 ) and summation over /i converges in the limit of e — > 0. 



Therefore, e^f becomes 



1 



2 (2tt)' 7o+ 
With use of the following formula 



roc r2kF 

/ / dq z ^2 

JQ+ J -2k F u 



Aire 



t2 



q\ + (q z + h) 



2kp 



exp 



I 2 



E 



sinhx 



n— — 00 



x 2 + (y + 2im) 2 2x (cosh x- cosy)' 

we can accomplish the summation over h as 

1 d 2 sinhd s (7x 
?i + (fts + ^) 2 2d s g ± (coshd s gj_ - cos d s q z ) 



E 



Consequently we obtain the following result: 



27r d s jo 



ds 



2tt/3 



df 1 



sinh s 



27r^ 7 cosh s — cos t 



exp 



s 

2d! 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



where d s = d s /ls denotes the inter-layer distance in units of the magnetic length. Figure ^ shows as a function 
of the inter-layer distance. We remark that e^f decreases with decreasing d s . This result originates from the fact 
that the inter-layer Coulomb interaction becomes more effective as layers get closer. 
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C. Comparison of Ground-State Energies 



We first compare the Coulomb energy per particle of the three states. Figure || shows the results. The energy of 
the (111) state is obtained by the VMC calculation for 5 layers with 42 particles in each layer. For the multi-Laughlin 
and the itinerant states the result is shown for the infinite-size system. We remark the energy of the multi-Laughlin 
state does not depend on the inter-layer distance since all layers are independent of each other and electrons feel other 
layers neutral. In the region of d s ^ Ib the Coulomb energy of the multi-Laughlin state is the lowest. This is evident 
because the system becomes a pile of almost independent two-dimensional layers in such a region. In the region of 
d s ^ Ib the (111) state has the lowest Coulomb energy. This result originates from the fact that the (111) state gains 
the inter- layer Coulomb energy at the cost of the intra-layer one in the region where intra- and inter-layer interactions 
are comparable. With respect to the itinerant state, the Coulomb energy decreases in the region of d s <, Ib- There is 
however no region where it is the lowest. 

So far we have considered only the Coulomb energy. We must take into account the tunneling energy of the itinerant 
state in order to compare the total energies of the three states. Let us assume that the transfer energy t depends on 
the inter-layer distance as 

t = Aexp (-ad s ) , (4-14) 

where A and a are constants. Then as shown in Fig. |?] we find a region of d s where the itinerant state has the lowest 
total energy. Here we put appropriate values into A and a, namely A = 0.2 e* 2 /ls, a — 1.0 /Ib for a case shown as 
ti, and A — 0.3 e* 2 /Ib, a = 0.6795 /Ib for another case shown as ti- The parameters for ti with B = 10 T are chosen 
to reproduce the band width of At = 2.5 meV = 0.179 e* 2 /Ib at d s = 0.226 A = 2.8 Ib which seems appropriate for a 
GaAs/AlGaAs superlattice 0. 

On the basis of the comparison we propose a phase diagram as shown in Fig. || for the ground state of the present 
system in the plane of t and d s . If we assume the relation ( 4. 14| ) between t and d s in the diagram, possible states 



follow a curve in Fig. || as d s is varied. Depending on the parameters A and a, there are either single transition fa) 
or double transitions fa)- Consequently the phase diagram Fig. [8] is interpreted as follows. 

(i) In the region d s ^> Ib, the system is a pile of independent two-dimensional layers and the incompressible 
multi-Laughlin state is stable. 

(ii) In the region of d s & Ib and the tunneling energy is small (t 0.1 e* 2 /ls), the (111) state which is also 
incompressible is stabilized instead of the multi-Laughlin state by the inter-layer correlation. 

(iii) When the tunneling energy is large, the compressible itinerant state is stabilized by competition of the tunneling 
energy and the Coulomb energy. 



V. DISCUSSION 



We compared the result of the (111) state for a finite size system and that of the multi-Laughlin and the itinerant 
states for the infinite size system. The finite size (N± and N z ) correction should be taken into account in the (111) 
state for more quantitative discussion. While we have no information with respect to the dependence on N z , the 
following remark is given from Fig. || about the dependence on N±. When one increases N± the Coulomb energy 
much increases in the region of d s 3> Ib- While for d s ~ Ib the Coulomb energy looks almost independent of iV_i_ 
within the error bars. Since the boundaries of the (111) state occurs with d s /ls < 2 in Fig. ^, the finite-size correction 
should not change the results qualitatively. 

It should be possible to observe transitions in actual systems by adding a pressure to the system or tilting the 
magnetic field. Since adding a pressure corresponds to control of d s , a change of parameters along a curve like t\ or £2 
in the diagram Fig. || is expected. On the other hand tilting the magnetic field corresponds to control of t independent 
of d s because the horizontal component of the magnetic field works on electrons to localize in the vertical direction 
and reduces t. Then the parameter change along a vertical straight line at a certain d s in Fig. || is expected. 

The transition between incompressible and compressible states is analogous to the Mott transition in the Hubbard 
model. Namely the incompressible multi-Laughlin and (111) states correspond to an insulator, while the compressible 
itinerant state corresponds to a metal. The transition between incompressible and compressible states occurs due to 
competition of the transfer (tunneling) and the on site (intra-layer) Coulomb energy. In our model we estimate that 
the critical transfer energy is almost 0.1 e* 2 /Ib- 

Although we discussed only three trial states for v = 1/3 in the present study, we cannot exclude the possibility 
of other states. In terms of two-dimensional incompressible states for v = 1/3, it is difficult to construct states 
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more stable than the multi-Laughlin and the (111) states within the simple family of Jastrow-type functions. One 
may however ask whether other states with three-body correlation or even more correlations are more favorable 
energetically. For three-dimensional compressible states, the itinerant state is taken to be the simplest homogeneous 
state. Then inhomogeneous states like CDW [^0|,^l) should also be discussed for more refined argument. Furthermore 
we have neglected the spin degree of freedom which may be relevant to some systems. For example the SDW state, 
which accounts for the QHE of (TMTSF) 2 X JIJ], can be a relevant state depending on parameters in the system. 
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APPENDIX: VMC METHOD FOR THE (111) STATE 

We describe the method of numerical calculation for the (111) state in this appendix. The expectation value of the 
Coulomb energy for the (111) state is represented by 



Ne (1 



J \ l p / \l=l' p<p' Kl' p,p< J 

-\ I (nn*^) ee v ^ - «tap}» i 2 i*^/3 m P m i 2 
v i p j i,i' p,p' 

= V!+V 2 . (Al) 

We can obtain the numerical value of £(m) for hnite size (finite particles) systems by the VMC method. In evaluating 
the second term with minimum statistical errors, we separate the variables in it and reduce the expectation value of 
the quantity between two bodies to that of one body. For this purpose we introduce the Fourier transform of the 
Coulomb potential V, one body distribution function /, and its Fourier transform g as follows: 

V (r - r>) = — / dnL% w (gj.) e "^( r ^), (A2) 



(27T) 



4-7re* 2 

Vzz- (q±) = cxp 



L„ . (\z — z 



(A3) 



/ {jlU d ^y^p- r ^\^i/sf' ( A4 ) 



f(r±) 



dill-) = / dr±f(rx)e~ 



/(nrw) e---ie i } 3 r, 



(A5) 



where we wrote z/d s as z. Then the second term V2 in (Al) is arranged into 



1 1 r°° 

V 2 = - -— 5 / Nl^22nq ± dqxVu> (qx)g(q±f 

Z ( 27r ) JO l y 



INN 1 ^ f°° ~ o 

oTTTT^E / ^q±dq ± V ll (q±)g(q±f. (A6) 
1 ( 2ir ) 1=1 J o 
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We use this formula to obtain the numerical value of V2 ■ Figure |9j shows an example of numerical results. The integral 
with respect to q± is performed by the trapezoidal rule because only discrete values of g (q±) are obtained by the 
VMC calculation. We cut off the upper region in the integral at a finite value q c because the integrand 

converges to zero near q± ~ q c and its contribution to the integral is negligible in the region q c < q±_ < 00. 
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FIG. 1. Configuration of the infinite-layer system. Layers are arranged perpendicular to the z-axis with the inter-layer 
distance d a . A magnetic field is applied along the z-direction. 

FIG. 2. Coulomb energy per particle of the multi-Laughlin state. The four points correspond to the particle number per layer 
N± = 20, 30, 42 and 72. The error bar indicates the Monte Carlo error. We estimate the value in the thermodynamic limit by 
extrapolating the linear line which fits the four points with the least squares method. The value is —0.4089 ± 0.0002 e* 2 /Ib- 
The energy does not depend on the inter-layer distance in our multi-Laughlin state. 

FIG. 3. Coulomb energy per particle of the (111) state vs. the inter-layer distance. The calculation is done for five layers 
with 20, 30 and 42 particles in each layer. 

FIG. 4. Common regions where two step functions, 9 (kF — \k\) and 6 (kF — \k — q z + g\), are nonzero. We have common 
regions in the cases (see text): (2) < q z — q < 2/cf and (3) —2kF < q z — g < 0. But we do not have a common region in the 
cases: (1) q z — g > 2&f and (4) q z — g < — 2fci?- 



FIG. 5. Coulomb energy per particle of the itinerant state in the thermodynamic limit. 
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FIG. 6. Comparison of Coulomb energies e per particle for the three states. The result on the (111) state is shown for a finite 
size system (N± = 42 and N z = 5), while the other two results are shown for the infinite size system (N± = oo and N z = oo). 
The (111) state has the lowest energy in the region of small inter-layer distance while the multi-Laughlin state becomes stable 
for large d s . The crossing of the two states occurs near d s ~ Ib- The itinerant state has no region where its Coulomb energy 
is the lowest. 

FIG. 7. Energy per particle of the itinerant state with the tunneling energy. We assume the dependence of t on the inter-layer 
distance as t = ^4exp (— ad s ) and show two exemplary cases. 

FIG. 8. Proposed phase diagram in the plane of t and d s . The solid lines are drawn on the basis of the VMC calculation 
for the (111) state of the finite size system (N± — 42 and N z = 5). The boundary with the itinerant state is determined from 
comparison with the Hartree-Fock calculation in the thermodynamic limit. 



FIG. 9. An example of g (q) calculated by VMC method. Calculation is performed for the system with N± = 20 and N z — 5. 
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